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Abstract 

One of the primary aims of the ongoing nuclear collisions at Relativistic Heavy Ion Collider (RHIC) and Large 
Hadron Collider (LHC) energies is to create a Quark Gluon Plasma (QGP). The heavy quarks constitutes a unique 
probe of the QGP properties. Both at RHIC and LHC energies a puzzling relation between the nuclear modihcation 
factor Raa{Pt) and the elliptic flow V 2 (pt) related to heavy quark has been observed which challenged all the existing 
models. 

We discuss how the temperature dependence of the heavy quark drag coefficient can address for a large part of such 
a puzzle. We have considered four different models to evaluate the temperature dependence of drag and diffusion 
coefficients propagating through a quark gluon plasma (QGP). All the four different models are set to reproduce the 
same Raa{Pt) experimentally observed at RHIC energy. We have found that for the same Raa{Pt) one can generate 
2-3 times more V 2 depending on the temperature dependence of the heavy quark drag coefficient. 


1. Introduction 

Heavy Quarks, charm and bottom, created in ultra- 
Relativistic Heavy Ion Collisions (uRHIC) represents 
ideal probes to study the Quark Gluon Plasma QGP 
III] El. An essential feature in analyzing Heavy quarks 
motion in a QGP is that their mass is much larger than 
the typical momentum exchanged with the plasma par¬ 
ticles entailing that many soft scatterings are necessary 
to change significantly the momentum and the trajec¬ 
tory of the heavy quarks. Therefore the propagation of 
heavy quarks has been usually treated as a Brownian 
motion that is described by means of the Fokker-Planck 
(FP) equation. In such an equation the interaction is en¬ 
coded in the drag and diffusion coefficient. The two key 
observables related to HQ that have been measured in 
experiments are the nuclear suppression factor Raa and 
the elliptic flow V2 Il3ll4l|5ll6jl71. Several theoretical 
efforts have been made using the Fokker Planck equa- 
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tion to reproduce the Raa and the V 2 experimentally ob¬ 
served ID El [ini [m Ell El El US Another approach 
used to describe heavy quark propagation is the Boltz¬ 
mann approach (BMl lfT6l[r7l[T8l[T9ll20ll2T]| . Analyzing 
the Raa gives information on the average magnitude of 
the interaction. We point out that studying the relation 
between Raa and V 2 it is possible to deduce other infor¬ 
mations on the interactions between HQ and bulk. All 
the approaches show some difficulties to describe simul¬ 
taneously Raa and V 2 - We have found that two ingredi¬ 
ents assume a particular importance in reducing the dif¬ 
ferences between experimental data for Raa and V 2 and 
theoretical calculations: the temperature dependence of 
the interaction and the use of full Boltzmann collision 
integral to study the time evolution of HQ momentum. 
The proceeding is organized as follows. In section 2 we 
discuss the Fokker Planck approach which is used to de¬ 
scribe the propagation of heavy quark through the QGP. 
In this section we present four different modelings to 
calculate the drag and diffusion coefficients. In section 
3 we compare the result for Raa and V 2 obtained with FP 
approach with those obtained using the BM approach. 
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2. Fokker-Planck approach 



Figure 1: Variation of the drag coefficients with respect to temperature 
at p = O.lGeV 


The propagation of the Heavy quark is described by 
the Fokker Planck equations which is indicated in the 
following equation 


dl 

dt 


dpi 


Mp)f + ^[Bijip)] 
dpj 


( 1 ) 


where A, and Bij are the drag and the diffusion coeffi¬ 
cients. To study the propagation of HQ this equation is 
replaced by the relativistic Langevin equation, which is 
more suited for numerical simulations 

dx: - —dt, 

E 

dpi = -Apidt + (V^) 

where t/x, and dpi are the coordinates and momenta 
changes in each time step dt, A is the drag force and 
Bq and B\ are respectively the longitudinal and the 
transverse diffusion coefficients; pj is the stochastic 
Gaussian distributed variable; Pfj - 6ij - PiPjlp^ and 

P/y - PiPjlP^ are the transverse and longitudinal pro¬ 
jector operators respectively. We employ the common 
assumption, Bq - Bi - D ll9l fTOl fT2l fT3l . To solve 
the Langevin equation a background medium which de¬ 
scribes the evolution of the bulk is necessary. In our 
work the evolution of the bulk is provided by a 3-i-lD 
relativistic transport code tuned at fixed rjjs Ealia 
which is able to reproduce the same results of hydro- 
dynamical simulations. The transport code provides at 
each time step the density profile and the temperature 
profile of the bulk which are necessary to calculate the 
drag coefficients. We have performed simulations of 
Am -h Am collisions ?A yfs - 200 AGeV for the mini¬ 
mum bias. In both cases the initial conditions in co¬ 
ordinate space are given by the Glauber model, while 



Figure 2: Comparison of the experimental data for the nuclear sup¬ 
pression factor as a function of pj at RHIC with the results we 
get with the FP equation for the four different T-dependences of the 
drag coefficient (Black lines). The orange lines show the results we 
get using the Boltzmann approach for the pQCD case and for the 
(^QPM{T),mq = = 0 case. 


in the momentum space a Boltzmann-Juttner distribu¬ 
tion function up to a transverse momentum pr GeV 
has been considered. At larger momenta mini-jet distri¬ 
butions, as calculated within pQCD at NLO order, have 
been employed Il24ll . The HQ distribution in momentum 
space is in accordance with the one in proton-proton that 
have been taken from ll25l . We have used four different 
modelings to calculate the drag coefficient entailing dif¬ 
ferent temperature dependence of the interaction . The 
diffusion coefficient is instead calculated in accordance 
with the Einstein relation D — YET. Our purpose is to 
investigate the relation between Raa and V 2 for differ¬ 
ent temperature dependence of the energy loss. In the 
first modeling we have evaluated the drag coefficient 
from (pQCD) and we have considered elastic interac¬ 
tion among HQ and the bulk (light quarks and gluons). 
The scattering matrix related to these processes MgHQ, 
MqHQ and AiqHQ in leading order are the well known 
Combridge matrix. The infrared singularity in the t- 
channel is regularized introducing a Debye screening 
mass iriD = yfAnal T with a running coupling ll26l . 

In another modeling we have evaluated the drag force 
from the gauge/string duality lIZTll through the follow¬ 
ing equation 


conf 


= c 


T2 

^ QCD 
Mr 


(3) 


where C = = 2.1 + 0.504)1. 

We have also considered another modeling in which the 
drag coefficient is evaluated considering a bulk con¬ 
sisting of particles with a T-dependent quasi-particle 
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Figure 3: Compaiison of the experimental data for the elliptic flow 
as a function of pj at RHIC with the result we get with the FP equa¬ 
tion for the four different T-dependences of the drag coefhcient (Black 
lines). The orange lines show the results we get using the Boltzmann 
approach for the pQCD case and for the Q!qpm{T), nig = mg = 0 case. 


masses, nig - ll3g^T^, nig = 3l4g^T^. This model 
is able to reproduce the thermodynamics of lattice 
QCD ll28l (see also ll29ll^[3n ') by htting the coupling 
g{T). Such a fit leads to the following coupling If28l : 


sHt) 


4871^ 

[(lliV,-2iV^)/n[T(^-^)]2 


(4) 


where A-2.6 and TITs-0.51. 

Finally we have considered a model in which the light 
quarks and gluons are massless but the coupling is from 
the QPM as indicated in Eq. This last case is indi¬ 
cated in the hgure as {uQPMiT), nig - trig - 0) and has 
to be considered as an expedient to have a drag which 
decreases with the temperature. 

For all the four cases considered the interaction has been 
rescaled to reproduce the Raa observed in experiments. 
We have simulated HQ propagation with the Langevin 
dynamics for the four different models presented above. 
The Langevin equation gives as output the momentum 
distributions of HQ at the quark-Hadron transition tem¬ 
perature Tc- The momenta distributions are convoluted 
with the Peterson fragmentation functions of the heavy 
quark indicated in Eq. [^in order to get the momentum 
distribution of D and B mesons. 


f(z) oc 


[z[l - i - I^]^] 


(5) 


where = 0.04 for charm quarks and = 0.005 for 
bottom quark. In hgure|^the nuclear modihcation factor 
Raa of the D and B mesons is shown as a function of pr 
for RHIC (200AGey). Instead in figure the elliptic 
flow (V2 = {{{pi - Ppiipl + Py))» at the same energy 
as a function of pj is depicted. We observe that the 


larger is the interaction in the region of low temperature 
the larger is the elliptic flow. The same conclusions has 
been discussed also in the light flavor sector as shown in 
Refs. ||32][33l. The reason of such a strong dependence 
of the elliptic flow on the temperature dependence of 
the drag coefficient is due to the fact that the elliptic 
flow is generated in the final stage of the evolution of 
the hreball when the temperature is lower. 


3. Boltzmann approach 

The Boltzmann equation for the HQ distribution 
function is indicated here 

P^d^fQ{x, p) = C[fQ\{x, p) ( 6 ) 

where C[/q](jc, p) is the relativistic Boltzmann-like col¬ 
lision integral which is solved by means of a stochastic 
algorithm. In such an algorithm whether a collision 
happen or not is sampled stochastically comparing the 
collision probability P 22 - VreiCTg^g+Q^g^g+Q ■ AtjAx 
with a random number extracted between 0 and 
1 OH EH EH [361. We use the Boltzmann equation to 
describe the propagation of the heavy quark and the 
evolution of the bulk 

The comparison between LV and BM approach has 
been studied in these references ClKIll where it is 
shown that for charm quark EP deviates significantly 
from the BM and such a deviation signihcantly depends 
on the the values of the Debye screening mass, whereas 
for bottom quarks the FP is a very good approximation. 
We considered in references nHiiii three values of 
nio: 0.4 GeV, 0.83 GeV and 1.6 GeV. Here we have 
not considered a fixed value of the Debye screening 
mass but a value which depends on the temperature 
according to nip, = gt. In hgures[^and[^the comparison 
for the Raa and V 2 at RHIC between the BM (orange 
lines) and the FP (black lines) are shown. We found 
that using the BM for the same values of the Raa we 
get larger values for V 2 with respect to those obtained 
using the FP. 

With the Boltzmann approach using the 
{<XQPM{T),mg - trig - 0) we get a value of the 
elliptic flow even larger with respect to the experimen¬ 
tal data, however this case represent an extreme and not 
realistic case. Our results show a non-negligible impact 
of the approximation in the collision integral involved 
in the Fokker Plack equation on the relation between 
Raa and V 2 - We summarize the results introducing 
a new plot in Fig. in which Raa vs V2 at a given 
momentum {pp - 1.3 GeV) is shown. This figures 
clearly shows how the building up of the V 2 can differ 
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Figure 4: Raa vs V2 obtained with the FP for the four different 
T-dependences of the drag coefficient with the experimental data at 
RHIC energy at pj = 1.3G^F. The open symbols indicate the results 
obtained using the BM approach 


up to a factor 3 for the same R^a depending on the 
temperature dependence of the interaction and on the 
approach, BM or FP, used to describe the propagation 
of the heavy quark in the QGR 


4. Conclusions 

We have found that different temperature depen¬ 
dences of the interaction can lead to difference in the 
elliptic flow V 2 by 2-3 times even if leading to the same 
Raa- Our studies suggest that the correct temperature 
dependence of the drag coefficient may not be larger 
power of T (as in pQCD or AdS/CFT) rather a lower 
power of T or may be constant in T. 

Moreover we have studied the difference in the build¬ 
ing up of the elliptic flow between Fokker-Planck and 
Boltzmann approach for a fixed Raa- We observe that 
the Boltzmann approach generates a larger elliptic flow 
than the Fokker-Planck lfT8l[T9l . 
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